One- and many-body effects on mirages in quantum corrals 
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Recent interesting experiments used scanning tunneling microscopy to study systems involving 
Kondo impurities in quantum corrals assembled on Cu or noble metal surfaces. The solution of 
the two-dimensional one-particle Schrodinger equation in a hard wall corral without impurity is 
useful to predict the conditions under which the Kondo effect can be projected to a remote location 
(the quantum mirage). To model a soft circular corral, we solve this equation under the potential 
WS{r — ro), where r is the distance to the center of the corral and ro its radius. We expand 
the Green's function of electron surface states for r < ro as a discrete sum of contributions 
from single poles at energies — iSi. The imaginary part Si is the half- width of the resonance 
produced by the soft confining potential, and turns out to be a simple increasing function of e^. In 
presence of an impurity, we solve the Anderson model at arbitrary temperatures using the resulting 
expression for Gg and perturbation theory up to second order in the Coulomb repulsion U. We 
calculate the resulting change in the differential conductance AdI /dV as a function of voltage and 
space, in circular and elliptical corrals, for different conditions, including those corresponding to 
recent experiments. The main features are reproduced. The role of the direct hybridization between 
impurity and bulk, the confinement potential, the size of the corral and temperature on the intensity 
of the mirage are analyzed. We also calculate spin-spin correlation functions. 



INTRODUCTION 

In recent years, the advances in scanning tunneling mi- 
croscopy (STM) made possible the manipulation of sin- 
gle atoms on top of a surface and the construction of 
quantum structures of arbitrary shape Ij. In particu- 
lar, quantum corrals have been assembled by depositing 
a close line of atoms or molecules on Cu or noble metal 
(111) surfaces IMIEQIS' These surfaces have the prop- 
erty that for small wave vectors parallel to the surface 
a parabolic band of two-dimensional (2D) surface states 
uncoupled to bulk states exists 6]. A circular corral of 
radius 71.3 A was constructed depositing 48 Fe atoms 
on a Cu(lll) surface 0]. The measured differential con- 
ductance dl/dV could be fitted by a combination of the 
density of eigenstates close to the Fermi energy ep, of a 
2D electron gas inside a hard wall circular corral ^2| . On 
another experiment, a Co atom acting as a magnetic im- 
purity was placed at one focus of an elliptical quantum 
corral, and the corresponding Kondo feature in dl/dV 
was observed not only at that position, but also with re- 
duced intensity at the other focus 4] . This "mirage " can 
be understood as the result of quantum interference in 
the way in which the Kondo effect is transmitted from 
one focus to the other by the different conduction eigen- 
states of a hard- wall ellipse ,7j. Among these eigenstates, 
the density of the one closest to is clearly displayed 
in the change of differential conductance AdI/dV after 
adding the impurity 0. More recently, the extrema in 
the degenerate S?"* and 38*'' conduction eigenstates of 
a hard wall circular corral have inspired an experiment 
with two impurities inside the corral |^. 

The theories of the mirage experiment can be classi- 
fied into those which start from eigenstates of a confined 



2D electron gas 0, 0, H E3, El E3 and those in which 
the confining atoms are modeled by a phenomenological 
scattering approach jl^ J^J . The latter have the disad- 
vantage that many-body effects are very hard to include. 
In addition, some features which are clear from the hard 
wall eigenstates, like mirages out of the foci 0,00,0], 
are somewhat hidden in the scattering approaches. On 
the other hand, while the eigenstates inside a hard wall 
corral are perfectly defined, in the actual experiments 
the boundaries of the corrals are soft and the eigenstates 
become resonances with finite width 6i 0|. This width 
plays a crucial role in the line sha pe o f AdI/dV and its 
magnitude at the mirage point f\ lllj. If 6i is large the 
mirage disappears, while if Si — and other parameters 
as in the experiment, there is no Kondo resonance ed ep- 
Ordinary Lanczos calculations have Si = and are 
unable to describe the line shape of AdI/dV. This can 
be corrected using an embedding method ^3- To the 
best of our knowledge, this theory and the perturbative 
one 0, ^1 using a constant value for Si = S as a parame- 
ter, are the only many-body ones presented so far, which 
can explain the voltage dependence of AdI/dV. Thus, a 
careful study of the broadening Si is necessary. 

Another important parameter which is still not well 
known is Sb/Sg, where St {Sg) is the part of the res- 
onant level width of the impurity state, which is due 
to hybridization with bulk (surface) states. Some the- 
ories have assumed that St ^ Ss, while others have taken 
'^fc = O000|T3,[iil. Clearly, S^ 0, because other- 
wise AdI /dV would be independent of the 2D conduction 
eigenstates in contrast to the experiment 0]. Instead, 
the line shape for both, the elliptical corral and the clean 
surface can be explained taking S^ ~ Q and the same set 
of parameters [ll|. A recent experiment suggests that 
Sb > Sg on the basis of the rapid decay in AdI /dV as the 
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STM tip is moved away from the impurity 0. How- 
ever, the 4s and 4p states of the impurity atom were not 
included in the analysis and they can affect the distance 
dependence of AdI/dV. Recently, it has been suggested 
that experiments with two impurities inside a quantum 
corral can elucidate the magnitude of Sb/Ss because the 
effect of interactions between impurities is expected to 
be proportional to S'^ ■ 

In this work we solve the impurity Anderson model 
which describes the mirage experiment using perturba- 
tion theory in the Coulomb repulsion U. The previ- 
ous theory of one of us J\ is extended to include non- 
vanishing db and a more realistic Si. We calculate the 
widths of the conduction states Si for a circular corral. 
We obtain an integral expression for the conduction elec- 
tron Green's function in the absence of the impurity, 
and then we expand it as a discrete sum of contributions 
from sirnple poles, using ideas borrowed from nuclear 
physics 01 and scattering theory 0|. This expansion 
takes a suitable form for the perturbative approach. We 
also discuss the space dependence of the compensation of 
the impurity spin. The calculation of is presented in 
section II. In section III we describe the model and the 
many-body approach. Section IV contains the results 
for AdI /dV in circular and elliptical corrals with a mag- 
netic impurity. In Section V, we calculate the correlation 
functions between the impurity spin and that of the con- 
duction electrons as a function of position. Section VI is 
a discussion. 



THE CLEAN CIRCULAR CORRAL 

Starting from free electrons in 2D, we model the 
boundary of an empty circular corral (without magnetic 
impurities inside) by a continuous potential WS{r — tq), 
where r is the distance to the center of the circle and tq 
its radius. The approximation of a continuous bound- 
ary (instead of discrete adatoms forming the corral) sim- 
plifies considerably the mathematics and is justified by 
the fact that the Fermi wave length for surface elec- 
trons 2TT/kp ~ 30 A is larger than the average dis- 
tance between adatoms ^ 10 A. The angular momen- 
tum projection perpendicular to the surface 1^ — m be- 
comes a good quantum number, due to the rotational 
symmetry around the center of the corral. For each m 
and energy E = (?ifc)^/(2m*), where m* is the effec- 
tive mass, the eigenstates can be written in the form 
'■Pkmi'f',9) = ipkmikr) e'™^ in polar coordinates, and the 
2D radial Schrodinger equation becomes 



r 19 



2mt r dr 



dipk 



dr 



+ [-^-^+WS{r-ro)-E]^krn = 0. 



(1) 

Its solution can be written in the form ipkm = '0fc~m^(^o ~ 
r) -|- ip^jy^d{r — ro), where 9{r) is the step function and 



(except for a normalization constant): 



^tm = A„,{k)J„,{kr)+B„,{k)Y„,{kr), (2) 

where Jm{Y„i) is the Bessel function of the first (second) 
kind. From the continuity of "iA/tm and the discontinuity in 
the first derivative according to Eq. , and using known 
expressions for the derivatives of the Bessel functions JJj , 
one obtains: 



Am{k) = 1 + 



2mlWl{h^k) 

Y-„+i(fcro) J„+i(/cro) ' 
y™(fcro) J™(fcro) 



Bm{k) = {l-A,n{k)) 



Yrn{kro)' 



(3) 



We normalize the wave functions (pkm {t, 0) inside a 
hard wall box of large radius R and take the limit 
R — > -|-oo. A quantity of interest in the many-body 
problem (see Eq. \ibl ) is the Green's function for surface 
conduction electrons in the absence of the impurity: 



km \ 1 J 

km 



Z ^kr\ 



(4) 



where the bar over i^km denotes complex conjugation and 
the sum over k extends over all radial wave vectors for 
which i^kmiR, d) — 0. Using Eqs. Q and ©, and taking 
the limit R +oo, we obtain after some algebra for 
r, r' < ro 



G2(z;r,0,r',0') = E 



dk k Jm{kr)Jm{kr') e™(«-»') 



2^(A2„(fc)+i?2Jfc))(, 



2m* 

(5) 

When W — > -l-oo, it is easy to see that Am,Bm oo, 
except for those values of k for which Ji{kro) — 0, and 
one recovers the known result for the hard wall problem 



J,„(7™r)J„(7™r') e™(«-«') 



(V7rroJm+i(7™ro))^ (2 



2m* 



(6) 



where 7™ro is the n*^ zero of Jm{ct)- For finite W one 
expects that the poles of Eq. (jBJ become resonances de- 
scribed by appropriate complex poles of the scattering 
matrix |l7l Il8l | . Actually, the integral of Eq. I^l can be 
thought as an integral over energy E{k) = (?lfc)^/(2m*) 
and for any z in the upper half of the complex plane, 
the integrand is analytical in E, except at the zeros of 
A"^ + . The physical zeros of A"^ + lie in the lower 
half plane of E. Then in principle one can evaluate the 
integral by the method of residues. A technical problem 
is that the Bessel functions diverge for infinite complex 
k and one has to assume a high energy cut off. This 
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FIG. 1: Partial contribution of conduction states with m = 
1 to the local density of states, and approximation using a 
summation of a finite number of contributions from poles in 
Eqs. Q, JHJ, as a function of energy for W = 15ft^/(2m*ro) 
and r = 0.14 ro. 



however has no physical consequences. Thus, can be 
written in the form: 



2 - e™ + i5"' 



(7) 



where fc™ are the complex zeros of A^{k) + B'^{k), with 
{hk-^ fl{2ml) = e;^-i(5™, 5™ > 0, and from the residues 
of the integrand in Eq. one obtains: 



a(^^(fc)+gi(fc)) I 

dk \''=K 



(8) 



In practice, in Eq.Q we include only the terms for which 
e™ < Ec, where the cut off energy Ec is typically taken 
as three times the Fermi energy. Thus, Eq. Q is an 
asymptotic low-energy expansion. It reduces to Eq. © 
for a hard-wall corral. 

A quantity of interest (because it is proportional to 
dl/dV lialSd) is local density of conduction states 
Ps{r, 9, e) which for r < ro is given by 



lim 



-lmG°{e + ir];r,e,r,e). (9) 

TT 



In Fig. ^we compare the contribution to the exact 
of the states with angular momentum projection m — I, 
with the corresponding expression using the partial sum- 
mation in Eq. 0. We have taken an effective mass 
m* = 0.38 me where me is the electron mass 0, l2lj . 
As in recent experiments i , we have taken tq = 63.5 
A, so that the fourth state (ordered in increasing en- 
ergy) with m = 1 falls ai ep Finally we took 
W = 15?i^/(2m*ro) ~ 2.37 eVA in order that the width 
of this state would be Si = 0.0213 eV, similar to that 
reported in some experiments 0- The point of observa- 
tion r = 0.14ro was taken at the maximum of Ji{jlr). 
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FIG. 2: Width of the energy levels as a function of energy, 
for (a) < C < 1 eV and (b) < e™ < 10 eV. 



One can see that the result in Fig^Jis consistent with the 
fact expected from previous studies with potentials that 
vanish except in a finite region [l7'| : the discrete sum Eq. 
(0) reproduces the low-energy part of G°. 

In Fig. 13 we represent (5™ as a function of e™. As a 
first approximation, (JJ^/e™ is a constant, independent of 
m and n. In any case, 5™ seems to depend on n and m 
only through e™. This result is surprising. For a given 
confining potential W, one would naively expect that 5™ 
increases with the radial part of the kinetic energy, being 
rather independent of the angular part. However, as m 
increases, the weight of the wave functions near the bor- 
der of the corral {r tq) also increases and compensates 
the decay in radial velocity. 

To study the effects of changes in the magnitude of 
the confining potential, we started form a situation with 
el = ep for rigid walls {W +oo). As W decreases, el 
decreases slightly, its width Si increases and the position 
rM of the maximum ofp2(r,0,e) for e = el deviates from 
its value 0.138ro for W +oo towards larger values. 
These changes are illustrated in Table HI The change 



in el can be interpreted, in a first approximation, as a 
small increase of the effective mass. In fact, using known 
properties of the Bessel functions 0|, it can be shown 
that for W ^ oo, the effective mass increases by a factor 
1 -I- h^/ (mlroW). Fitting the lowest energies of a circular 
corral with a model with hard walls, gives m* = 0.38 TOg 
I23. For moderate W, it seems that Si ~ l/W, how- 
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W 



Sl [meV] e\ [meV\ Maximum [% To] 



0.1 


171.4 


364.8 


14.81 


1 


100.2 


381.3 


14.74 


5 


51.0 


393.1 


14.60 


7 


41.2 


396.1 


14.55 


10 


31.3 


400.7 


14.48 


12 


26.6 


403.1 


14.44 


15 


21.3 


406.4 


14.39 


18 


17.4 


409.3 


14.34 


oo 





441.1 


13.80 



TABLE I: Variation of the energy e\, its width 5\ and the 
position of the maximum of p{r,6,e\) as a function of the 
confining potential W . 



ever for large W, i5™ ~ l/W"^ . In contrast to Sl, is 
only weakly dependent on W (see Tabled). 

What happens with e™ and (5™ if the size of the corral 
is changed? We expect that W remains constant, inde- 
pendent of tq. If To is increased by a factor /, changing 
variables in the Schrodinger equation JQ), to p = r/ f, 
one realizes that (except for normalization) the following 
equality holds for the wave functions of given projection, 
position, energy and parameters of the potential: 



V>mir, e, e, fro, W) = ^m{j,0. ft, tq, fW). 



(10) 



Then, increasing by a factor / one expects, as a first 
approximation, a radial dilation by a factor / of the cor- 
responding wave function, a compression of the energy 
levels e by a factor of f , and for large W , that the widths 
5i remain approximately constant. 



THE MANY-BODY PROBLEM 

In this section, we explain the model and approxima- 
tions used to describe the electronic structure of a system 
composed of a quantum corral and a magnetic impurity 
inside it. Basically we consider surface states, described 
inside the corral by a Green function like that dis- 
cussed in the previous section, hybridized with one lo- 
calized d orbital with an important on-site repulsion U 
0, 0, 0, 0| ■ We also include an hybridization of the d 
orbital with bulk states. We are neglecting the degener- 
acy of the d level, which we believe is not important for 
the essential physics, and other s,p orbitals brought by 
the impurity which might affect the magnitude of dl /dV 
near the site impurity. 

The Hamiltonian can be written as: 



H 



JO' 



EdY^dl 



+uY,d\d^dld^ + VAyj{R,)dUj, + H.c] -f 
+ E^b[4V + H.c.]. (11) 



where sj^. {bjcr) £^re creation operators for an electron in 
the j*'' surface (bulk) conduction cigenstatc in the ab- 
sence of the impurity but including the corral. The im- 
purity is placed at the two-dimensional position Ri on 
the surface, and we assume that the hybridization of the 
impurity d orbital with the surface state j is proportional 
to its normalized wave function at that point (pj{Ri). We 
write the proportionality constant as VgX, where Vs is an 
energy (representing the hybridization in a tight binding 
model [9|,lll|) and A = 2.38 A is the square root of the 
surface per Cu atom of a Cu(lll) surface. 

At sufficiently low temperature, the differential con- 
ductance dl/dV when the STM tip is at position Rt is 
proportional to the density of the mixed state (2C| . 



fcriRt) = XYipjiRt)Sja +qda 



(12) 



dm 

dV 



PficF + eV) 



— ImG^eF + eV), (13) 



where G^(a;) — {{fa',fa))uj is the Green function of 
fa{Rt), and q{Rt) is the ratio between the tunnehng ma- 
trix elements tip-impurity and tip-surface, and is signifi- 
cant only for very small Rt — Ri, due to the rapid decay 
of the wave functions for d states. 
Using the equations of motion 



(cj-ej)((s 
(w-ej,)((d^; 



jo; dl^))u 



- S,,, + VsXip^{R,)({dy,s],J) 
= VAiPj(R,){{dy,dl))^, 

- F,A^„(i?.)((rfo;4»-, 



(14) 



is expressed in terms of the Green function for the 

d electrons G'^{uj) = ((do-; and the unperturbed 

conduction electron Green function for surface states G° 
. We drop the superscript a in the following because of 
the spin independence of the problem. The difference 
in G/ between the results with and without impurity 
becomes: 



AG/ = {Vs\^Gl{Rt,Rr,uj) + q) X 
X {VA^G'i{R„Rt,u:) + q) Gd,. 



(15) 



where 



GO(i.„E„.)=E^^^^^^^^- (16) 



JO 



We assume that G^ is known from the one-body prob- 
lem. For the circular corral, it has been discussed in the 
previous section. 
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The remaining task is to calculate Gd using a many- 
body approach. The most accurate calculation would 
be to use the Wilson renormalization group How- 
ever, the particular structure of the one-body problem, 
and the lack of continuous symmetries in the general 
case renders its application very difficult. We have used 
perturbation theory up to second order in U Eif . 
Near the symmetric case Ed + U/2 = e^?, the theory is 
quantitatively correct up to U ^ 8A, where A is the 
resonant level width [22|. Out of the symmetric case, 
interpolative self-consistent schemes Esl Eoj still work 
for moderately large values of U js^, l3l| . In particular, 
the persistent current in small rings with quantum dots 
practically coincides with exact results [IJl- Here we as- 
sume a situation in which the effective one-particle d level 
E'^^^ = Ed + U{dl^dcr) is very near the Fermi level. This 
is justified by first-principle calculations fs^l and allows 
us to avoid selfconsistency. 

Then: 



[G,H]-i = [G°H]-i-SH 



(17) 



where G° is G 



for U — and Ed replaced by E'^J^ ^ and 
is the self-energy up to second order in IJ evaluated 
from a Feynmann diagram involving the analytical ex- 
tension of the time ordered G\ to Matsubara frequencies 

[Ulii: 

I](zc^,,T) = V''TY,Gl{iu:i-ivM^Vra), (18) 

XilVra) = -r^G°(iw„)G°(ic^„+il/„), (19) 

n 

where LOn = (2rt + X^ixT and = Itl-kT . 

From the equation of motion of Gd for ?7 = and using 
Eqs. lfn|l we obtain 



where we have replaced 



UJ + irj — ei 



, (20) 



(21) 



assuming that the V^^ and the bulk density of states are 
featureless near ep and can be replaced by a constant. 
The real part of the sum can be absorbed in E'^-^-^ . 

From the discussion of section „ and previous works 
01, we know that G^{Ri, Ri,uj) can be expanded as a 
sum of contributions from discrete poles. In particular, 
for the circular corral, G^{uj;r,d,r' ,9') is given by Eqs. 
Q) and JSJ. For other shapes of the corral, the results of 
the previous section suggest that within errors of a few 
percent: 



Gg (w, Ri, R2) 



E 



^'=^{R,)r,{R2) 

ui — Cj + idpij/^F^ 



(22) 



where now (p'j{R) are the discrete eigenstates of the hard 
wall corral and ej are their energies, calculated with a 
slightly renormalized mass, dp is the width of the reso- 
nance at the Fermi level. In practice we cut off the sum in 
Eq. H22|l at an energy ~ Sef, retaining N poles. Then, 
G^{u!) can be expressed as a sum of + 1 poles and 
residues at these poles [3^ Denoting the poles of G"((jj) 
by ri — iSi and their residues by a;, we can write: 



N+l 

E 

1=1 



ai 



bj -ri + i6i 



(23) 



This is the retarded Green function. In order to evaluate 
perturbative diagrams, hke Eqs. ((TH|l . l(T^ . one needs the 
analytical extension of the time ordered Green function 
to imaginary frequencies lu —> iun |3 Issf : 



G°(zc.„) = J2 



ai +ai + {ai - Q;)sgn(a;„) 
2 [iujn + iSiSgn{uJn) - u] 



(24) 



where TiJ is the complex conjugate of a/ and sgn(a::) is 
the sign of x. Here and in what follows, the origin of 
energies has been placed at the Fermi energy (ep = 0) 
for simplicity. 

Using standard methods 'si'], the sums over Matsub- 
ara frequencies, Eqs. H18|) . H19|) can be transformed into 
integrals over branch cuts. The first integral can be done 
analytically in terms of the digamma function ^-(2) 
A lengthy but straightforward algebra after analytical 
continuation back to real lu leads to: 



S(w) 



where 



1 



1 



GS(c.-y)Imx(2/) 



1 



ev/T + 1 



imG^d{~v) x(^ + y)], 



(25) 



Ztt ^ — ' 



Im 



- {BimH + Airni-iv)) 4';(w) - 

-(b;;;{-uj) + a^,{lu))Wi{-lo), (26) 



with 



ai a,i 



UJ + ri - Trn - iSl + i6n 

LU ~ n + rrn + i6i + i6„ 



T ^ ^ T /'I + *(n — ^) 



(27) 



Summarizing, our approximation for the Green function 
of the d electrons Gd consist in the following steps: we 
first decompose the unperturbed Green function G^ given 
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by Eq. and or into a sum of poles, Eq. (j^ . 
Then, the self-energy is calculated using Eqs. 

and ^7^, and finally is obtained from Eq. ((T7jl . 
The change in dl/dV after adding the impurity can then 
be evaluated using Eqs. l(T^ and l(T^ . 



RESULTS 
The elliptical corral 

The line shape of the change in differential conductance 
after adding the impurity, AdI/dV has been reported for 
the elliptical corral with eccentricity e = 0.5 and semi- 
major axis of 71.3 A with a Co impurity placed at the 
left focus 0. Comparison with our results determines 
the value of Vg which leads to the observed line width 
~ 0.009 eV. We took U = 1 eV and assumed a situa- 
tion near the symmetric case, for which the Hartree-Fock 
E^/^ level is near ep 113 ■ We shifted E"/^ so that the 
minimum of AdI/dV as a function of energy lies near the 
experimental position. E^^-^ might be calculated using 
self-consistent interpolative approaches f [2^ EgHsollsij ) . 
but we avoided them here. We have taken two values for 
the level width of the conduction states at ep, Sp — 0.02 
eV, suggested by some experiments |3| and 6p = 0.04 
eV, which leads to the observed ratio 1/8 between the 
intensity in dl/dV at the mirage {—Ri) and at the im- 
purity {Ri) For 5p = 0.04 eV and 6b = 0, we obtain 
that Vs = 0.68 eV and q ^ 0.03, leads approximately 
to the observed line shape. They also explain the line 
shape when the impurity is placed on a clean surface 
[m . We then decrease the value of Vs by a factor of \/2 
to Vs = 0.48 eV and increase the resonant level width 
of the impurity to Sb = 0.032 eV in such a way that the 
width of the Kondo peak in the impurity density of states 
Pd{^) (not shown) is the same as before. This conditions 
implies approximately the same contribution from bulk 
and surface states to the Kondo resonance, as assumed 
in previous studies |l3l |. As shown in Fig. |31 and in 
agreement with previous calculations both sets of 
parameters lead to very similar results except for a scale 
factor. This factor is approximately Vg , as expected from 
Eq. (0 

The decrease in the intensity of the resonance at the 
mirage point with respect to the impurity position can 
be understood from Eqs. H15|l and H22|l and is a conse- 
quence of the destructive interference of the 42*'* state 
of the ellipse (which lies at the Fermi level and is even 
under reflection through the minor axis a) with other 
states which are odd under a (like states 32, 35, and 51): 
all terms contribute to the same sign to the dominant 
imaginary part of the sum in Eq. (|22(l for Ri — i?2, 
while there is a partial cancellation for i?i = — i?2- This 
interference decreases with decreasing 5p and the in- 
tensity at the mirage point increases. For 5p = 0.02 
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FIG. 3; Change in the spectral density of the / state, pro- 
portional to AdI/dV (see Eqs. 11311 . 11511 ) as a function of 
frequency at the impurity site for two values of q and at the 
other focus, for two sets of parameters (a) 14 = 0.68 eV, 
St = 0, E'/^ = 0.015 eV (b) K = 0.48 eV, 5^ = 0.032 eV, 



V 

T = 0. 



Ey^ = 0.007 eV. In both cases Sp = 0.04 eV, f/ = 1 eV, 



eV, Vs should be decreased shghtly to 0.63 eV to keep 
the same width of the resonance, for Sb = 0. Taking 
Vs — 0.45 eV and increasing Sb from zero to 6b — 0.035 
eV to keep the width of pd{(^) (as above), the result for 
A / dujpf{uj){-df{uj)/du) - AdI/dV is shown in Fig. H 
for two different temperatures. At T = 0, in comparison 
with the previous results for Sp — 0.04 eV, the inten- 
sity at the mirage increased to nearly 60% of that at the 
impurity position. At the Kondo temperature, the reso- 
nances at both positions broaden and the amplitude at 
the minima decrease more than 50% of the zero temper- 
ature result. Most of the broadening is due to the effect 
of the derivative of the Fermi function, —df{Lu)/duj. In 
fact, Apf{uj) at the impurity looses only ^ 1/3 of the 
intensity and broadens very little as T is increased from 
to 0.005 eV. 

If 6p is further reduced, the space dependence of 
Apf(uj) is almost exactly that of the conduction eigen- 
state at ep and therefore the intensities at both foci are 
practically the same. However, the shape of the reso- 
nance strongly changes, and two positive contributions to 
Apf{Lu) appear at both sides of ei;-. This has been shown 
before for Sb = 0, 0| and is a consequence of the splitting 
of the Kondo peak in the impurity density of states pd{u>) 
0, 01 ■ This splitting of the Kondo peak into two peaks 
at both sides of the Fermi energy, which was obtained 
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FIG. 4; Change in differential conductance at both foci for 
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FIG. 6: Impurity spectral density for several values of 5f and 
T indicated in meV inside the figure. Other parameters as in 
Fig. H 
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FIG. 5: Apfcr ~ AdI /dV as a function of frequency for Sf = 
0.001 eV and other parameters as in Fig. 2] 

first in perturbation theory I?, Ill and later by exact cal- 
culation of a reduced Hamiltonian plus embedding in the 
rest of the system ,12j has been recently confirmed by 
numerical Wilson renormalization group calculations in 
a simpler system with U(l) symmetry As shown in 
Fig. |S1 this effect remains for 6h 0, although the split 
peaks in Pd{^) do not reduce to two delta functions for 
(5f = in this case (see Fig. O. In Fig. Owe also show 
the temperature dependence of Pdi^) for the parameters 
of Fig. El 

So far, the experiments were done at very small tem- 
peratures (T ^ 4 K ^ Tk — 50 K), and there is no 
appreciable difference with the T — results. A study 
of the temperature dependence would confirm that the 
resonance is due to a many-body effect. In general (in- 
cluding other situations discussed below), the scale of the 
T dependence is determined by Tk, which in turn is given 
by half of the with of the peak in the impurity spectral 
density Pdiy)- As can be seen comparing Figs. 4 and 6, 
for moderate Sp (hke 5f = 0.02 eV), this width is larger 



than the width of the dip in the conduction electron spec- 
tral density pdi^) or in pf{ijj). 

The essential features of the space dependence of 
AdI / dV are very similar to that discussed previously 
0, 01 and will not be reproduced here. Essentially it 
reflects the electron density of the conduction state at 
ep, somewhat blurred out of the impurity position for 
6f > 0.03 eV. 

In the following we use the four set of parameters 
which lead to a reasonable line shape at the impurity 
site {Sf = 0.02 eV or 6f = 0.04 eV, 6b = or 6b ^ 0) to 
study other situations for which the voltage dependence 
of AdI/dV has not been reported. We first consider an 
elliptical corral with the same eccentricity but semima- 
jor axis reduced from a = 71.3 A to a = 64.6 A. In 
this case, the conduction state 35 (instead of 42) of the 
ellipse falls at e^. The electronic density of this state 
1^35(^)1^1 has maxima at i? = ±0.4a, displaced from 
the foci R — ±0.5a. Studying the negative interference 
at the mirage point, one expects a stronger mirage at 
Rt — 0.4a if the impurity is placed at Ri = —0.4a 
We studied two cases, with the impurity at one focus or 
at one maximum of |(p35(i?)|^ for the four set of parame- 
ters mentioned above (actually, from the discussion of the 
previous section, one would expect a 10% reduction of 6f, 
with respect to the larger ellipse, but we neglect it). For 
Ri — —0.5a (impurity on focus), the resonance narrows 
with respect to the previous case. This is mainly due to 
the fact that the effective hybridization of the state at 
the Fermi level is proportional to |(p35(i?i)|^ which is ap- 
proximately 60% of \(p42{Ri)\^ of the larger ellipse. The 
total width is near 0.005 eV for 6b = , and near 0.007 
eV for (5b ^ (see Fig. H (a)). The intensity of the de- 
pression at the other focus Rt = 0.5a is roughly half of 
that at Rt = 0.4a at the extremum of ip35{Ri) near that 
focus. In particular for 6f = 0.02 eV (not shown), the 
dip in Apf^Lu) at Rt = 0.4a is slightly more pronounced 
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FIG. 7: Ap/CT ~ AdI/dV as a function of frequency at several 
points for an ellipse with a — 64.6 Aand two impurity posi- 
tions: (a) Ri = —0.5a and (b) Ri = —0.4a. Other parameters 
as in Fig 2] 



than that at the impurity position Rt = Ri = —0.5a. 

When the impurity is placed at one extremum of 
'PasiR)^ Rt = 0.4a, the peak broadens with respect to 
the larger eUipse to - 0.015 eV for db = 0, or 0.012 eV 
for Sb ^ 0, due to the ~ 20% larger effective hybridization 
with the state at e^? (see Fig. 13(b)). As expected 7|> 
the intensity at the mirage point Rt — 0.4a, compared to 
the impurity point is larger in the smaller ellipse. 

While the intensity of the mirage is closely related with 
the hybridization of the impurity with the state at e^?, 
the width of the peak in the impurity spectral density and 
hence Tk depends on the effective hybridization with the 
impurity of all conduction states. For example, if in the 
ellipse of the original experiments one suppresses artifi- 
cially the hybridization with the main odd states 32, 35, 
and 51, the width in Pd{^) is reduced to roughly half its 
value for realistic Sp, while the magnitude of /S.dI/dV 
at the mirage point increases substantially, reaching al- 
most the same magnitude as in the focus occupied by the 
impurity. The width in AdI/dV is less sensitive to the 
parameters as that of Pd{'-^)- As expected, this sensitivity 
is decreased for increasing 5p. 



The circular corral 

In this section we study the voltage and space depen- 
dence of AdI/dV for a circular corral, for conditions cor- 
responding to recent experiments 0: a radius rp ~ 63.5 
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FIG. 8: Ap/o- ~ AdI /dV as a function of frequency for Rt = 
Ri = 0.146ro (dashed line) and Rt — — O.lSro (solid line) for 
the circular corral. The thinner lines are the results using 
the approximate expression (1221 changing Ri to 0.138 ro and 
Rt = ~Ri. Parameters are W = 7n^ / {2mtro) , V = 0.48 eV, 
St, = 0.032 eV, E^/-'' = 0.01 eV and q = 0.02. 



A, such that the degenerate 37*'' and 38*'' conduction 
eigenstates lie at the Fermi energy. In contrast to the el- 
liptical corral, we do not use the approximate expression 
l)22|l for the conduction electron Green's function, but 
the correct expansion given by Eqs. and ©. The 
width of the resonant states is determined by the poten- 
tial scattering at the boundary W, and its value at the 
Fermi energy Sp is given in Table 

In Fig ISl we show Ap/o- ~ AdI/dV for an impurity 
placed at a distance 0.146 ro to the center (correspond- 
ing to a maximum of the conduction density of states, 
see Table m, for W — 7?i^/(2TO*ro), and other parame- 
ters similar to Figs. EJb) or [3 As in previous cases, an 
intense resonance is observed not only at Ri but also near 
—Ri. The maximum of the absolute value of Apf„ for the 
mirage is slightly displaced towards O.lSro. We also show 
in Fig |S1 the results using the approximation for the un- 
perturbed conduction electron Green's function based 
on the wave functions for the hard wall, Eq. 1)22(1 . The 
resonance is broader in this case. This is due to the fact 
that with decreasing W , not only the maximum in real 
space of the conduction density of states shifts, but also 
the spectral weight of the resonance at near this max- 
imum decreases, leading to a smaller effective hybridiza- 
tion. When W is increased so that Sp ^ 0.02 eV, the 
intensity at —Ri increases. Also when Jf, = is taken, 
increasing Vs the resonances broaden by ~ 20%. These 
changes are qualitatively similar to those reported for the 
smaller ellipse. 

The space dependence of the conduction density of 
states is shown in Fig. The two clear minima near 
±i?i can be seen, one of them corresponding to the im- 
purity and the other to the mirage, in agreement with 
experiments . The addition of the impurity breaks the 
degeneracy of the states with m = ±1 at e^- If the zero 
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FIG. 9: Conduction electron density of states at ep as a 
function of position. Other parameters as in Fig. |H1 



line of the angular variable 9 is aligned with the impurity, 
the wave function which hybridizes with the impurity has 
an angular dependence proportional to cos 9, and the re- 
maining one, proportional to sin 9 does not hybridize. 
This has led to the proposal 5j that two independent and 
simultaneous mirages can be observed placing two impu- 
rities at the same distance from the center and forming an 
angle of 90 degrees. While this seems in agreement with 
experiment, the conduction states with even m (which 
lie out of e^), introduce a small interaction between both 
channels. 

To end this section we discuss what happens when rg 
is slightly reduced in such a way that the degenerate 35*'' 
and 36*'' conduction eigenstates lie at the Fermi energy. 
These states have m = ±3 so that one expects five dif- 
ferent mirages as a function of the angle 9, located near 
r = 0.32ro, where the eigenstates for a hard wall corral 
are peaked. Since the value of the wave function at these 
extrema is nearly 3/4 times smaller than in the previous 
case, the effective hybridization of the impurity with the 
states at the Fermi energy is reduced. As a consequence, 
the Kondo depression observed in the energy dependence 
of AdI/dV is narrower (with a total width of ~ 0.008 eV 
for ~ 0.02 eV and Si, ^ 0), and the intensity at the 
mirages is smaller. In Fig. 1101 we show the space de- 
pendence of Apsa ^ /S.dl / dV for parameters near those 
of Fig. 4. A potential W corresponding to harder walls 
(leading to Sp ^ 0.02 eV) than in Fig. 9 has been cho- 
sen in order that the different mirages can be clearly seen. 




FIG. 10: Contour plot of the conduction electron density 
of states at EF as a function of position for a smaller corral 
(see text). Parameters are W = 157i^/(2m*ro), V = 0.45 eV, 
& = 0.035 eV, and E^" = 0.013 eV. 

The most intense ones are those for 6* = ±120 degrees. In 
addition to the five depressions in Apsa{Rt) correspond- 
ing to the mirages, one can see a ring of positive Ap^o- 
centered at the impurity site at Ri = 0.337ro and radius 
smaller than tq. The intensity of this ring increases with 
decreasing W (softer walls) and is therefore also appar- 
ent in Fig. O We believe that this feature is probably 
related with the Friedel oscillations which are observed 
experimentally when the impurity is placed on the clean 
surface [l5| . 



SPIN-SPIN CORRELATIONS 

The conventional view of the Kondo effect for an im- 
purity in an ordinary metal, interprets it in terms of a 
magnetic screening cloud around the impurity of radius 

C = Tivf/Tk, (28) 

where fp- is the Fermi velocity. The existence of this 
cloud is still controversial Recent theoretical 

work has shown that the persistent current as a func- 
tion of flux j($) in mesoscopic rings with quantum dots 
changes its shape smoothly as the length of the ring L 
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FIG. 11: Function g{r /£) (see Eq. Il34ll 'l for different symmet- 
ric limits of integration indicated inside the figure. 



goes th rough ^ and that jL is a universal function of i/^ 

In spite of this controversy, it is interesting to study the 
effect of the corral on the space dependence of spin-spin 
correlations. To this end we consider the correlations 
between the impurity spin S and that of the surface con- 
duction states at position r, s(r) 



(S-s(r)) 



E 



(29) 

For infinite confinement {W +oo, dp = 0), exact nu- 
merical results in a reduced basis set show that the space 
dependence of this function follows closely the probabil- 
ity of the wave function lying at the Fermi level [l^ . 
However, as discussed before, this situation is not realis- 
tic. 

In the uncorrelated case U = 0, using Wick's theorem 
and the symmetry of the ground state one has: 



(S-s(r))=--A 



(30) 



The saddle point approximation of the slave boson treat- 
ment of the Anderson model , or of the 1 /N expan- 
sion of the Kondo model js^ , lead to a similar expression. 
Here we use Eq. (|30ll as an approximation for the finite 
U problem. This amounts to an infinite partial summa- 
tion of diagrams which can be separated into two pieces. 
The expectation value entering Eq. H30|l can be related 
to Gd{w) without further approximations using equations 
of motion. At zero temperature this leads to: 



(S-s(r))-- 



3A^ 
27r2 



3 



dwlm- 



Gdito) 



IT] - 

(31) 



We believe that the main qualitative features of the prob- 
lem are retained by this approximation, at least for r ^ ^ 
which is the relevant case for the corrals. To check this, 
we evaluated this expression for the case of a Kondo im- 
purity in an ordinary three-dimensional metal, assuming 
(as usual in this case) that plane waves describe the con- 
duction electrons, and (except for a multiplicative con- 
stant) Grf(w) ~ l/(w-eF + irK) 00, with Tk < tF- 
Calling litj) the integral entering Eq. H31(l and lineariz- 
ing the spectrum around tp = hvphp one has: 

'^lpj{Ki)ipj{r)I{ej) ^ J (fkexp{ik ■ r)I{hvFk) 



kp f 

— / dk sm( kr) I (hvpk). 
r J 



(32) 

Here the vectors are denoted by boldface (k, r) to distin- 
guish them from their absolute values (fc,r). Changing 
variable z — (fc— kp)^, and performing the integral in lo 
{I{hvpk)), we obtain except for a constant: 



, sin^ (kpr) , , , , , 

(S.s(r)) i^5(r/f), (33) 



g{r/C) = J dzcos(— ) ^2 _^ / ■ (34) 

g{0) diverges unless one introduces a cutoff at \z\ = 
kpS, = cf/Tk- For other values of the argument the 
integral can be extended to infinity at both sides, and it 
can be related to tabulated integrals The result is: 



9{x) 



2^Ei(-a;); x 7^ 0, 



(35) 



where Ei(x) = J -00 '^^ 6xp(i) /t is the integral exponential 
function. The function g(r/^) is represented in Fig llll 
for several values of the cutoff in |z|. As seen in the 
figure, this cutoff only modifies the function for small 
arguments. There is a crossover for r ~ ^ and for r ^ 
g{r/^) decays linearly with r, leading to a decay 
of (S • s(r)). The functional form of (S • s(r)) is very 
similar to that of the spin susceptibility suggested by 
Sorensen and Affleck |37ll38l | for r < ^. However, for r ^ 
^ they obtain an exponential decay. Therefore, we expect 
that the main features of the space dependence of the 
correlation function are retained by the approximation 
Eq. (1201) for r < ^. 

In Fig. 1121 we show the evaluation of Eq. (|31() for the 
elliptical corral, with the parameters of Fig0] and zero 
temperature. In contrast to dl/dV or the local spectral 
density of states near ep (Fig. IH and previous calcula- 
tions for perfect confinement {5p = 0) [13, the spin cor- 
relations do not recover substantially at the empty focus. 



11 




FIG. 12: Correlations between the spin of the impurity and 
that of the conduction electrons as a function of the position 
of the latter. Other parameters as in Fig. |1] 

although still faint features related with the conduction 
state at the Fermi level (j = 42) are still discernible.. 
The ratio of the spin correlations at both foci is ~ 1/170. 
Note that from Tk and vp for surface electrons in the 
mirage experiments, one gets ^ 90 A, slightly larger 
than the distance between foci. As it is clear from the 
calculation above, for r < ^ all conduction states (also 
those at high energies) contribute to the spin-spin corre- 
lations. Therefore, the destructive interference near the 
empty focus is larger than in the case of dl / dV , in which 
mainly unperturbed conduction states near ep are sam- 
pled. The participation of high energy electrons in the 
formation of the singlet ground state has been stressed 
recently [s^. 

SUMMARY AND DISCUSSION 

We have studied the 2D Green's function for free elec- 
trons, subject to a potential that simulates a continuous 
circular corral. We disregarded the hybridization with 
bulk states at the boundary of the corral. In scattering 
theories, an imaginary part is introduced for the phase 
shift to take into account this effect in phenomenological 
way 0, 0, 0| , and this leads to loss of intensity at the 
mirage point when an impurity is introduced in the cor- 
ral. In our theory this loss of intensity comes as a result 
of quantum interference, as explained in Section IV. Al- 
though the density of states is continuous in space and 
energy, the energy dependence of the Green's function at 
low energies can be written as a sum of discrete poles, 
representing resonant states. This property is expected 
to persist for more realistic potentials which might in- 
clude hybridization with bulk states or individual atoms 



at the boundary at the corral [iJl . The introduction of a 
width Si of each resonance is the main difference with a 
calculation assuming a hard wall corral. The energies 
and the space dependence of the density of states inside 
corral are well described by the calculation assuming hard 
walls, but the magnitude of the latter is somewhat over- 
estimated. Our model leads to a width of each resonance 
5i, which is to a first approximation, proportional to its 
energy Ci . We must warn however that at least above the 
Fermi energy, the surface states have an intrinsic width, 
even in absence of the corral |^. Thus, while we expect 
that 5i can be reduced for example by depositing a sec- 
ond line of atoms at the boundary of the corral, there 
are probably intrinsic limits to this reduction. Another 
shortcoming of our one-body calculation is that actually 
the effective mass depends on energy. Photoemission ex- 
periments suggest a flattening of the dispersion above e p 
6]. This is also apparent in the comparison of experi- 
mental energy levels for a circular corral an those of a 
hard wall calculation Q and persists in our results for 
soft boundaries. 

The expansion as a sum of discrete poles of the con- 
duction electron Green's function is a convenient starting 
point for a perturbative treatment of the Anderson model 
that describes the physics when a magnetic impurity is 
introduced inside the corral. The line shape of the change 
in differential conductance AdI/dV and the relative in- 
tensity at points distant from the impurity (where a "mi- 
rage" of the impurity is observed j^) is very sensitive to 
Si. This sensitivity is reduced but persists if an impor- 
tant intrinsic width Sb of the impurity level due to direct 
hybridization with bulk states is introduced. A smaller Si 
leads to a stronger intensity at the mirage point, and the 
space dependence tends to that of the conduction state 
at ep- However, for very small Si the line shape (volt- 
age dependence), is strongly distorted and two peaks of 
positive weight of AdI /dV appears at moderate non-zero 
voltages. The line shape and its width depend also on 
the particular conditions of the experiment. This fact 
includes the case of the clean surface, as observed exper- 
imentally 0, and calculated before 0. A larger width 
and more intense resonance is expected when the wave 
function of the conduction state which lies at the Fermi 
energy has a stronger amplitude at the impurity. 

While the space dependence of AdI/dV is determined 
only by the conduction electron Green's function in ab- 
sence of the impurity (see Eqs. (|13|l and (|15(l '). the voltage 
dependence is very sensitive to the impurity Green's func- 
tion. Most previous theories either assume or fit this line 
shape dllElli or are unable to explain the observed one 
because Si = was assumed 0. Also while numerical 
diagonalization with Si = and a reduced basis set show 
spin-spin correlations which remind the conduction state 
at the Fermi level we obtain using additional ap- 
proximations discussed in Section V, that no appreciable 
projection to the mirage point is present in these corre- 
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lations. This is due to the fact that the screening of the 
spin involves conduction states far from the Fermi energy 
ep, while AdI/dV selects states near ep. These experi- 
ments seem unable to settle the still controversial issue 
of the Kondo cloud [s^jIUjII^- Instead, a crossover as a 
function of size is expected in mesoscopic rings |^ 0| . 

Our many-body theory is able to explain the main fea- 
tures of the space and voltage dependence of AdI/dV in 
several experiments involving the quantum mirage 0, Q 
and predictions for other geometries and experimental 
conditions were made. While it would be interesting to 
test this predictions and study the effects of temperature, 
the theory has some limitations. A real Co atom added 
to a Cu(lll) surface has degenerate localized 3d levels 
with strong correlations. We have taken only one d or- 
bital with a moderate on-site Coulomb repulsion U = 1 
eV to be able to use the perturbative approach. We be- 
lieve however that these shortcomings affect quantitative 
details but do not invalidate our conclusions. Previous 
results using larger values of U and a different technique 
led to results very similar to the perturbative ones |l2l| . 
A Co atom also has extended 4s and 4p electrons, which 
should have a large hopping to the STM tip when it is 
near the impurity. This should affect the amplitude of 
AdI/dV at moderate distances to the impurity. This ef- 
fect has so far been neglected in the theories of the mirage 
effect. 
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